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Abstract 

We have calculated the heat current in the normal metal/insulator/ supercon- 
ductor contacts with arbitrary transparency of the insulator barrier. In the 
tunneling limit (small transparencies), the heat flow out of the normal metal 
reaches its maximum at temperature T ~ 0.3 A. At higher values of trans- 
parency, the interplay between single-particle tunneling and Andreev reflec- 
tion determines optimum transparency which maximizes the density of heat 
flow out of the normal metal. In clean contacts, the optimum transparency is 
about 0.1 at T = 0.3A and decreases with temperature roughly as (T/A) 3 / 2 . 
In disordered contacts, disorder enhances Andreev reflection and shifts the 
optimum point towards smaller transparencies. The optimal ratio of the bar- 
rier resistance to the resistance of the normal electrode is Rn/Rt — 0.01 at 
T = 0.3A and decreases with temperature similarly to clean contacts. For 
disordered contacts we also plot current-voltage characteristics for arbitrary 
values of the ratio Rn/Rt- 
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I. INTRODUCTION 



The flow of electric current in the normal metal/insulator/superconductor (NIS) contacts 
is accompanied by the heat transfer from the normal metal into the superconductor. This 
principle can be applied to the refrigeration of electrons in the normal metal Jl|. (Implic- 
itly, the same principle is used in the enhancement of the superconductivity in the SIS'IS 
structures - see and references therein.) The mechanism of the heat transfer in the 
NIS contacts is the same as that of the well-known Peltier effect in metal/semiconductor 
contacts - see, e.g., ||. Due to the energy gap in the superconductor, electrons with higher 
energies (above the gap) are removed from the normal metal more effectively than those 
with lower energies. This makes the electron energy distribution sharper, thus decreasing 
the effective temperature of electron gas in the normal metal. 

In the limit of very small transparencies D of the insulator barrier, when the mechanism 
of electron transport across the barrier is the single-particle tunneling, the magnitude of heat 
flow out of the normal metal increases with transparency. However, at larger transparencies, 
coherent two-electron tunneling ( "Andreev reflection" ) starts to dominate electron transport 
and suppress the heat flow. This occurs because in the Andreev reflection electrons with all 
energies, including those inside the energy gap, are removed from the normal metal. Below 
we study the crossover between the two regimes and find the optimum transparency which 
maximizes the heat flow through a unit area of the NS interface. The calculations are carried 
out in the two cases of contacts with clean and disordered electrodes. It is shown that, in 
accordance with the general understanding that disorder enhances Andreev reflection , 
in disordered contacts the optimum transparency is shifted towards smaller transparencies 
with increasing disorder. 
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II. CLEAN CONTACTS 



The model of NIS contact we consider is a constriction between normal metal and super- 
conductor with characteristic dimensions d that are much smaller than the coherence length 
£ and inelastic scattering length in the electrodes. Because of the condition d -C £, we can 
neglect variation of the superconducting order parameter A in the vicinity of the constric- 
tion and solve the problem, assuming that A is constant in space up to the NS boundary 
and is equal to its equilibrium value inside the superconductor. 

The properties of such a constriction depend strongly on the relation between d and 
elastic scattering length I near the junction. In the clean limit {d <C £), electron motion in 
the constriction is naturally decomposed into several uncoupled transverse modes (provided 
that the NS interface is smooth and conserves transverse momentum). In this case one 
can follow Blonder, Klapwijk, and Tinkham || and solve the Bogolyubov-de Gennes (BdG) 
equations independently for each transverse mode. The basic result of such a solution (for 
details see, e.g., the review ||) is contained in the probabilities of normal (B) and Andreev 
(A) reflection from the NS interface as functions of the quasiparticle energy e: 

, n _ D 2 W)\* fl(l-2Rea 2 (e) + |a(e)| 4 ) 

U 1 - 2RRea 2 (e) + R 2 \a(e)\ 4 ' U 1 - 2i?Rea 2 (e) + i? 2 |a(e)| 4 ' {> 

Here D, R are transmission and reflection probabilities of the insulator barrier, D + R = 1, 

which are assumed to be independent of energy on the energy scale given by A, and a(e) is 

the amplitude of Andreev reflection from the ideal NS interface with D = 1: 



1 e - sign(e)(e 2 - A 2 ) 1/2 , | e |> A, 
= -r < (2) 
I e _,(A 2 -e 2 ) 1 /2 | e |< A . 



Using the reflection probabilities ([I]), we can write the balance equation for energy distri- 
bution of electrons moving to and from the NS interface |J. Energy distribution of electrons 
f + (e) that move from the bulk of the normal metal to the NS interface is the equilibrium 
Fermi distribution shifted by eV, / + (e) = /(e — eV). Electrons moving from the interface 
into the normal metal are produced in three processes: quasiparticles incident from the 



superconductor are transmitted into the normal metal with the probability (1 — A — B); 
electrons are reflected from the interface with probability B; and holes are Andreev-reflected 
as electrons with probability A. (The latter process can be described in more direct terms as 
tunneling of Cooper pairs out of the superconductor.) Thus, the energy distribution / - (e) 
of electrons moving into the normal metal is: 

/-(e) = A(e)(l - /+(-e)) + B(e)f+ (e) + (1 - A(e) - B(e))f(e) . (3) 

From the thermodynamic relation dU = dQ + /idN applied to electron gas in the normal 
metal the heat current j flowing from the normal metal into the superconductor in one 
(spin-degenerate) transverse mode can be calculated as 

3 = ^-Jde(e-eV)(f+(e)-f-(e)). (4) 

It is straightforward to check that when deviations of electron energy distribution from 
equilibrium in each electrode are small, eq. (f|) is equivalent to another frequently used 
expression for the heat current, j = Th, where T is the temperature of the electrode and h 
is the entropy flow. 

The last step in the calculation is a summation over transverse modes. We do this 
summation assuming that the electron motion in the junction is quasiclassical, and that NS 
interface is smooth and conserves momentum along it. Then the transmission probability D 
of the interface depends only on the energy e± of electron motion across the interface, and can 
be taken to be D(e±) = D Q exp{(e± — /i)/e tr }. In this expression, D is the transparency at 
Fermi energy, and e^ r is an energy scale associated with electron motion under the barrier, 
which is assumed to be much larger than the superconducting energy gap. Under these 
assumptions summation over transverse modes (i.e., over angles of incidence on the interface) 
gives for the total heat current J: 

iV r D o dD r 

J= ^hL ~D J d <t-z V )U + {t)-r^D)), (5) 

where N = S^netr/nTi 2 is the effective number of transverse modes, m is electron mass, and 
S is the junction area. 
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Taking the limit D — > in eqs. (P and one can see that for D — > eq. (|) for the 
heat current reduces to the form given by the tunnel Hamiltonian approach: 

3 - ik SZ * e(f2 - A2 >^fef^ £ - eV > - ^ • < 6 > 

where Rt is the normal-state tunneling resistance of the barrier, R^ 1 = Ne 2 Do/%h. 

Figure 1 shows the heat current as a function of the bias voltage V across the contact for 
different transparencies D calculated numerically from eqs. (||), (|3|), and (H). The curves 
illustrate how increasing transparency of the barrier suppresses heat transfer out of the 
normal metal. We see that for each transparency there is an optimal bias voltage which 
maximizes the heat current J. The inset in Fig. 1 shows J at the optimal bias voltage in 
the tunneling limit calculated from eq. (0). The heat current J in this limit is maximum 
at T ~ 0.3A and decreases both at small and large temperatures. One can check that at 
T< A the heat current decreases as (T/A) 3 / 2 . 

In Figure 2 we plot the heat current per one transverse mode at the optimum bias voltage 
as a function of barrier transparency Dq. At small transparencies the heat current increases 
linearly with transparency, indicating that we are in the tunneling limit where electron 
transport is dominated by the single-particle tunneling. However, at larger transparencies 
the heat current starts to decrease with transparency due to increasing contribution to 
transport from the two-particle tunneling. At the transition point between the two regimes 
the heat flow out of the normal metal is maximized. A crude estimate of the transparency 
which corresponds to this transition at low temperatures can be obtained if we compare the 
amount of heat per transverse mode generated by the two-particle tunneling in the normal 
electrode, ~ (DqA) 2 /h, with the heat flow out of this electrode in the tunneling regime. As 
was noted above, the latter can be estimated as (T 3 A) 1//2 Z) /fi (see also inset in Fig. 1). 
From these estimates we see that at T < A the optimal transparency scales as (T/A) 3//2 . 
This conclusion is in qualitative agreement with Fig. 2. 
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III. DISORDERED CONTACTS 



In this section we consider the same model of a short NIS constriction, but assume 
that the constriction dimensions d are much larger than the elastic mean free path I in 
the electrodes. A convenient way to describe such a disordered constriction is provided by 
the quasiclassical kinetic equations for non-equilibrium Green's function G of the electrodes 
|T0|^T2[. Green's function G is a triangular matrix in the Keldysh space: 

( qR qK \ 



G 



\ ° a j 



(7) 



where the advanced G A , retarded G R , and Keldysh component G K are 2x2 matrices in the 
electron-hole (Nambu) space. The retarded and advanced functions carry information about 
the excitation spectrum of the system, while the Keldysh function describes the distribution 
of quasiparticles. 

For small constrictions, d <C £, we can neglect all but the gradient term in the equation 
for Green's function in the vicinity of the constriction, so that the equation is reduced to 



the diffusion equation flQHT2 



V{VGVG) = , (8) 

with the diffusion coefficient T> = ^Vp£, where vf is the Fermi velocity. Diffusion equation 
in the electrodes should be complemented by the boundary condition at the NS interface 
H0: 



cnGWG = [6, G s ] ■ (9) 

In this expression, n is the vector normal to the NS interface, and G, G$ are the Green's 
functions on the normal and superconducting sides of the interface, respectively. The coef- 
ficient c describes the interface transparency and can be written as 2a /g, where o = 2e 2 z/D 
is the conductivity of the normal metal, and g = e 2 v(nvFD j 1 R) is the normal-state con- 
ductance per unit area of the interface. Here D and R are the transmission and reflection 
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probabilities of the interface, v is the density of states at the Fermi level, and (...) denotes 
averaging over angles of incidence on the interface. 

On a quantitative level, properties of the constriction depend on its specific geometry. 
Below we consider a simple one-dimensional (ID) model shown schematically in Fig. 3. In 
this model the constriction is represented as a ID normal conductor of length d, cross-section 
area A and resistance Rn = d/(Aa), and it is assumed that the contribution of the bulk 
regions to the total resistance of the structure is much smaller than Rjy. In this case we can 
neglect variations of the order parameter A and Green's function Gs in the superconducting 
electrode, and assume that A and Gs are equal to their equilibrium values up to the NS 
boundary: 



G 



R(A) 



' g R(A) fR(A) ^ 

K -f R{A) -g R(A) j 



1 



e ± iO) 2 - A 2 



-A -e 



(10) 



G K S = G R (e)n(e) - n(e)G A (e) , n(e) = tanh( 



2T' 



(11) 



In the ID constriction (Fig. 3) the Green's functions depend only on the coordinate x along 
the constriction, so that the diffusion equation becomes: 



dx ^ dx ^ 



(12) 



Normalization condition on G, G 2 = 1, [ |T2"| , |H)| implies that the retarded and advanced 
components can be parameterized as follows [TB|: 



G R(A \e) = ±(a z cosh[U R{A \e)] + ia y smh[U R(A) (e)]) 



(13) 



where a's are Pauli matrices, and upper and lower signs are for G R and G A respectively. 
The fact that the retarded and advanced functions are nondiagonal indicates the proximity 
effect in the normal region. U R( - A ^ (e) determines the local density of states as a function of 



the energy e. The parameterization (|13|) allows us to think about it as an imaginary polar 
angle in the 2-dimensional space of (a z , a y ). The magnitude of the angle corresponds to the 
degree of "proximization" of the normal region. 



The Green's functions should satisfy the equilibrium boundary condition at the normal 
end of the constriction (x = —d): G R ( A > = ±cr z , i.e., U R ^ = 0. Equation flT2"|) with this 
boundary condition determines U R ^ >{x): 

UR(A) =a R(A) {1 + X^ (M) 

where a R ( A \e) = U R( - A \x = 0,e). Substituting eqs. ([13|) and flU]) into the boundary condi- 
tion (0) we reduce it to a transcendental equation: 

± ^ a R ^ = f R ^ cosha*^ - g R ^ sinha^ . (15) 
Rn 

Here g R( ~ A ^ and f R ( A ) are components of the equilibrium Green's functions ([10|) of the super- 
conductor, and Rt is the normal-state tunneling resistance of the NS interface, R^ 1 = gA. 
The next step is to find the Keldysh function G K which can be represented as [ [T2|j : 



G = G R (e)h(e) - h(e)G A (e) , h(e) = A(e)l + f z (e)a z . (16) 

Equation ([12]) with the retarded and advanced components flllf ) gives the following equations 
for the distribution functions f^ z ) ([T6|) : 

(1 + cosh(C/* T f/A ))^ fl = const = B Kz)/d. (17) 

Combining these equations with the boundary condition ([J) at x = and equilibrium 
boundary condition at x — —d: /i,( z )(— rf, e) = (n(e — eV) ± n(e + eV))/2, we find the 
functions which determine currents across the NS interface: 

Bi(e) = ^(2n(c) - n(e - eV) - n(e + eV)) , B z (e) = ^-(n(e - eV) - n(e + eV)) , 

(18) 

where 

A 1{z) = (g R - g A )(cosha R + cosher 4 ) - {f R t f A )(sinha R ± sinha A ) , 

r Rt Mnh[(a R T a A )/2] 



Equations flT5|) and (|T8|) enable us to find electric current and heat current in the NS 
junction. The electric current can be written as follows: 

I = —— / de tr{a z [GVGf(e)} = — / deB z (e) . (19) 

Figure 4 shows the current / ( |iTS| ) as a function of the bias voltage V at several ratios of 
the resistances Rt/Rn and vanishing temperature. At large tunnel resistances Rt the I — V 
characteristic exhibits the gap at V < A/e due to vanishing density of states in the super- 
conductor, and associated singularity at V = A/e. At smaller Rt the gap is closed by the 
increasing contribution to the current from Andreev reflection, but the I — V characteristics 
still has the singularity (logarithmic divergence of the differential conductance) at V = A/e. 
Transition between the two regimes is clearly visible in the zero-bias linear conductance 
shown in the inset in Fig. 4. At Rt < Rn the conductance with a good accuracy equals 
simply (Rt + Rn)' 1 , while at Rt > Rn it decreases as Rn /R\. 

Recalling the definition of the heat current used in the previous section and eq. ( |19D for 
the electric current, we can write the following expression for the heat flow out of the normal 
electrode in disordered NIS junction: 

1 z'+oo 

J = j ^ de(e - eF)(S 1 (e) + B z {e)) . (20) 

As can be seen from eqs. ( |i~5|) and (|18D , the currents Bi z have the property Bi^(— e) = 
^Bu z \(e). With these relations, eq. fl2~0|) can be rewritten as follows: 

1 r°° 

J =- IV + -n^ deeB^e). (21) 
e z RN Jo 

At Rt 3> Rn, both eq. ([0J) for the electric current and eqs. (|20"D, (pT|) for the heat current 
are reduced to the form given by the tunnel Hamiltonian approach. In particular, the heat 
current J is again given by eq. (||). 

Some results of the numerical calculation of the heat current from eqs. (0), flHf), and 
(PT|) are shown in Figs. 5 and 6. We see that the properties of the heat current in disordered 
contacts are qualitatively similar to those for clean contacts. The heat current density grows 
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with increasing transparency of the tunneling barrier in the regime of small transparencies 
(large tunnel resistances Rt) and is gradually suppressed at large transparencies. The main 
difference with the clean case is that the scale of transparencies (in particular, the optimum 
transparency) is shifted downwards by a small factor £/d, which decreases with increasing 
resistance of the normal electrode. The physical reason for this shift is that disorder in 
the normal electrode enhances contribution of Andreev reflection to transport and therefore 
suppresses the heat flow out of the normal metal. 

IV. CONCLUSION 

In conclusion, we have calculated the heat current in clean and disordered NIS micro- 
contacts which is caused by electric current flow across the NS interface. The mechanism of 
the heat transfer is analogous to that of the Peltier effect in normal metal/semiconductor 
structures. Results for clean NIS junctions are obtained by solving the Bogolyubov-de 
Gennes equations. Disordered junctions are described with the quasiclassical equations for 
non-equilibrium Green's functions of the electrodes. In both cases the heat current density 
exhibits non-monotonic dependence on the interface transparency, increasing at small trans- 
parencies and decreasing at large transparencies. The transition between the two regimes 
takes place at the transparency which is determined by the interplay of single-particle tun- 
neling and Andreev reflection. 

At intermediate temperatures, T ~ A, this transition occurs at the transparencies which 
are larger than the barrier transparencies of the typical tunnel junctions between good met- 
als. For instance, even quite small specific tunneling resistance on the order of 10 Ohmx /im 2 
corresponds to barrier transparency ~ 10~ 4 . However, as it was demonstrated above, as tem- 
perature decreases, the transition point moves rapidly towards smaller transparencies, and 
at low temperatures finite barrier transparency can pose an important limitation on the 
refrigeration power of NIS junctions. As follows from our estimates in Sec. 2, in the clean 
case transparency-related limitation should become important at T ~ 0.01A in junctions 
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with specific resistance 10 Ohmx/im 2 or less. 

The authors gratefully acknowledge discussions with K. Likharev. This work was sup- 
ported in part by ONR grant # N00014-95-1-0762. 
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FIGURES 



Figure 1. Heat current J in the clean NIS contact versus bias voltage V for several 
transparencies Do of the insulator barrier calculated from eq. (^). From top to bottom, 
D = (tunneling limit), 0.03, 0.1, and 0.2. The inset shows the heat current in the 
tunneling limit calculated for the optimum bias voltage as a function of temperature; at low 
temperatures J oc (T/A) 3 / 2 . 

Figure 2. The maximum heat current density in the clean NIS contact as a function of 
transparency Dq of the insulator barrier for several temperatures. For discussion see text. 

Figure 3. Schematic diagram of disordered NIS contact. Darker region shows the quasi- 
1D constriction of length d which determines the resistance Rn of the normal electrode. 
An insulator barrier with normal-state resistance Rt is placed at the NS interface in the 
constriction. 

Figure 4. I — V characteristics of disordered NIS contact at zero temperature and several 
values of the resistance ratio Rt/Rn- From top to bottom, R T /R N = 0, 1, 3, 10, 100. The 
inset shows the linear conductance G = dl/dV \v=o as a function of Rt/Rn- The curves 
illustrate the transition between the tunneling regime (Rt > Rn) characterized by the gap 
at V < A/e and "metallic" regime (Rt < Rn) characterized by large subgap conductance 
and excess current at V ^> A/e. 

Figure 5. Heat current J in the disordered NIS contact versus bias voltage V for several 
ratios of the normal electrode resistance Rn to the resistance Rt of the insulator barrier. 
As in the clean contacts, the heat current is maximized at V ~ A/e. 
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Figure 6. The heat current density in the disordered NIS contact calculated at the 
optimum bias voltage as a function of the insulator barrier resistance Rt- The curves are 
similar to those in Fig. 2 for clean contacts. Note, however, that the scale of the x-axis, 
Rn/Rt, for disordered contacts corresponds to much smaller transparencies of the tunnel 
barrier than for clean contacts (R N /R T — (d/£)D 3> D ). 
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